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Abstract 


The  present  report  is  the  third  and  final  deliverable  accor¬ 
ding  to  EOARD  Contract  SPC-94-4081,  item  0003. 

The  paper  contains  the  dynamical  problem  solutions  for  certain 
thermal  modes  of  artificial  plasma  expansion  in  space,  codes  to 
calculate  the  parameters  of  freely  expanding  plasma  jet  and  plasmo- 
id,  representative  samples  of  calculations,  comparison  of  theoreti¬ 
cal  results  obtained  with  data  of  measurements  in  plumes  of  elect¬ 
ric  and  conventional  gas  thrusters,  comparison  of  self-similar  so¬ 
lutions  obtained  with  published  data  of  analytical  and  numerical 
calculation  of  parameters  for  plasma  and  gas  jets. 


Nomenclature 


a,  1  =  transverse  and  longitudinal  spatial  scales  of  plasma  jet 

e  =  Ion  charge 

k=(a’0)-2  =  factor  of  flow  divergency 

M  =  Mach  number 

n  =  dimensionless  criterion 

m  =  Ion  mass 

N  =  particles  flow  rate 

n  =  particles  concentration 

T  =  temperature 

x,  r,  <p  =  cylindrical  coordinates 

u,  v  =  x-  and  r-  components  of  flow  velocity 

V  =  flow  velocity 

Y  =  specific  heat  ratio 

X  =  mean  free  path 

6  =  plasma  electric  conductivity 

-ri  =  r/a  =  self-similar  variable 

v  =  numerical  multiplier 

C  =  separation  constant 

Subscripts: 

"c "  =  quantity  at  flow  axis 

"m"  =  maximal  quantity 

"0"  =  initial  quantity 

Prime  "<M  denotes  spatial  along  axis  x  derivatives. 
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Introduction 


The  objectives  of  the  work  is  development  of  dynamical  and 
mathematical  models  of  exhausted  plasma  plumes  and  plasmoids  in 
space  with  due  regard  for  actions  of  inner  electric  fields  and 
thermal  conductivity  in  plasma. 

These  models  are  intended  to  predict  the  behavior  of  artifici¬ 
al  plasma  formations  (APF)  arising  in  space  particularly  in  vicini¬ 
ty  of  a  spacecraft  as  the  typical  result  of  normal  operation  of 
certain  onboard  systems  including  electric  thrusters  both  steady 
state  and  pulse  operation  ,  plasma  contactors  for  tether  systems 
and  for  spacecraft  charging  control  as  well  as  during  active  expe¬ 
riments.  Knowledge  of  the  APF  characteristics  is  important  for  de¬ 
termination  of  communication,  navigation  and  ranging  conditions  , 
for  evaluation  of  APF  influence  on  spacecraft  environment  and  con¬ 
sequently  on  sensitive  onboard  equipment  operation  as  well  as  for 
active  experiments  planning  and  data  processing. 

The  previous  intermediate  report  submitted  as  the  second  deli¬ 
verable  according  to  EOARD  Contract  SPC-94-4081,  item  0002,  conta¬ 
ins  dynamical  problem  formulation,  general  and  particular  solutions 
describing  evolution  of  artificial  plasma  formations  injected  from 
a  spacecraft  and  a  brief  description  of  algorithms  and  code  under 
development  for  calculation  of  spatial  and  temporal  plasma  proper¬ 
ties  distributions  in  the  frame  of  Development  of  Dynamical  and 
Mathematical  Models  of  Exhausted  Plasma  Plumes  and  Plasmoids  in 
Space. 

The  dynamical  problem  describing  evolution  of  both  fully  ioni¬ 
zed  stationary  supersonic  plasma  jet  and  plasmoid  is  formulated  and 
solved  in  a  self-similar  form.  The  basic  idea  of  this  approach  is 
that  a  non-uniform  plasma  or  gas  flow  tends  "to  forget  its  histo¬ 
ry",  i.e.  initial  conditions,  so  that  the  latter  becomes  actually 
insufficient  for  flow  expansion  far  from  the  source.  Because  of 
subsequent  reasonable  simplification  the  input  combined  partial 
differential  equations  of  two  components  plasma  dynamics  were  rear¬ 
ranged  in  a  set  of  separable  equations  so  that  the  problem  was  re¬ 
duced  to  quadratures  and  to  solution  of  ordinary  nonlinear  diffe¬ 
rential  equations. 

For  practically  interesting  cases  of  plasma  jet  with  uniform 
temperature  or  enthalpy  in  the  flow  cross  sectional  plane  the  solu¬ 
tion  was  obtained  in  a  simple  analytic  form. 

The  most  noticeable  feature  of  three  dimensional  plasmoid  ex¬ 
pansion  is  the  inversion  of  its  shape. 


1.  Representative  thermal  modes  of  plasma  Free  Expansion 

1.1  Adiabatic  plasma  or  gas  jet  expansion 
In  the  previous  intermediate  report  [1]  the  dynamical  problem 
describing  expansion  of  fully  ionized  stationary  plasma  jet  is  for¬ 
mulated  and  solved  in  a  self-similar  form.  For  two-dimensional  su¬ 
personic  plasma  jet  of  low  divergency,  the  variables  x  and  r  in  the 
final  set  of  dynamical  equations  separates. 

The  self-similar  form  of  solution  for  this  problem  is  written 
as: 

u  =  uc  (x)  - y (n) :  v  =  u-da/dx-ri;  T  =  Tc  (x)  •  x (n) ;  nu  =  f(t\)-Nv  /tta2. 
where  r  =  a(x)  is  a  certain  flow  current  line. 

It  was  found  that  all  dissipative  effects  mentioned  above  may 
be  evaluated  using  the  only  dimensionless  criterion 
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Here  "a"  is  a  spatial  scale  of  plasma  jet  transverse  size. 

For  the  condition  of 
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the  combined  equations  (6)  of  report  [1]  may  be  rearranged  to: 
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Here  C  is  free  separation  constant  for  variables  x  and  n  =  r/a. 


The  equations  1.1,  1.2,  1.5  and  1.6  result  from  equations  of 

motion  for  coordinates  x  and  r,  expression  1.3  represents  equation 
of  energy  and  expression  1.4  is  the  separation  condition  for  vari¬ 
ables  x  and  i\.  The  equations  1.1,  1.3  and  1.4  yield  relationships: 

T/T0  =  1  -  (u2/u02  -  1) -M2  •  (1-1 )/2 
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Therefore  the  separation  condition  1.4  is  consistent  with 
equations  1.1  and  1.3  if 
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The  solution  obtained  is  valid  for  conditions  u  *  const  with 
M02  »  1,  i.e.  for  more  wide  range  of  jet  parameters  than  that 
presented  in  report  [1], 

Finally  the  following  Jet  parameters  distributions  for  M02  »  1 
are  obtained: 
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As  a/a0  »  1  or  for  M02  -►  oo  this  solution  is  identical  to  that 
presented  in  report  [1]  for  conically  expanding  Jet. 


1.2  Isothermal  expansion  of  ion  beam 

Operation  of  onboard  Stationary  plasma  (SPT)  or  Anode  layer 
(TAL)  electric  thrusters  results  in  arising  in  space  of  expanding 
ion  beam  with  the  following  representative  initial  parameters 
[3+6]  N  ~  2-  1019c-1 ,  T  -  1+5  eV,  6  =  103+104  1/Om-m,  a~0.1,  a‘~0.2 
corresponding  to  criterion  n  =  20  +  100. 

Laboratory  measurements  (see  e.g.  Ref. [5])  suggest  that  the 
electron  temperature  in  the  beam  is  practically  constant  ranging 
within  1+10  eV. 

For  this  isothermal  modes  of  plasma  free  expansion,  the 
solution  of  the  combined  equations  (6)  of  report  [1]  may  be  written 
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The  separation  condition  2.4  is  consistent  with  equation  of 
motion  2.1  if  the  following  relationships  are  satisfied: 
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Therefore  the  separation  condition  for  variables  x  and  i\  is 
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Finally  the  following  beam  parameters  distributions  are 
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The  initial  part  of  SPT  or  TAL  plume  is  characterized  (see  e.g. 
[4]  +[7]  )  with  the  following  conditions: 

')(M02=  muo2/T0  ~  500-100  »  1,  so  that  u/u0  =  1  and  a'=a0'. 
Resulting  parameters  distributions  in  conically  expanding  beam 
are  the  following: 

-  transverse  size:  a2  =  x2/k; 

-  factor  of  flow  divergency:  k  =  (a‘0)~2; 

ft  k 

-  concentration:  n  = - :  2.8 

2Jtuc  x2(l+r2/a§) 

eN  k 

-  beam  current  density  along  axis  x:  ox  =  enu  = - —  : 

2%  x2  (l+kr2/a2  )3/2 

mftuc  k 

-  impuls  flux  density  along  axis  x:  mnu2  - - : 

4Jt  x2  (l+kr2/a2  )2 

Transverse  parameters  distributions  in  isothermally  expanding 
beam  are  more  wide  than  that  of  adiabatically  expanding  jet  so  that 
peripherial  energy  amd  substance  fluxes  in  the  beam  are  relatively 
greater. 


2.  Program  and  Calculation  of  Plasma  free  Expansion 

2.  i.  General  Description  of  Program 

The  dynamical  problem  solution,  presented  in  report  [1],  is  a 
mathematical  base  for  development  of  PC  program  for  calculation  of 
plasma  free  expansion  in  space. 

PC  program  DYNAMIC  developed  consists  of  two  units  JET  and  PLAS- 
MOID  for  calculation  of  flow  characteristics  of  steady  state  and  pulse 
plasma  injections  correspondingly. 

The  former  calculates  and  plots: 

-  profile  of  plasma  jet  characteristic  transverse  size  along  the 
jet  axis; 

-  radial  dependencies  of  plasma  concentration,  temperature  and 
velocity; 

-  lines  of  equal  plasma  concentration. 

Input  data  are  the  particles  flow  rate  and  initial  values  of  tem¬ 
perature,  velocity,  characteristic  jet  size  and  its  derivative  given 
at  a  plasma  source  exit. 

Unit  PLASMOID  generally  calculates  three-  dimensional  plasmoid 
parameters  for  a  time  moment  fixed  resulting  in: 

-  time  profiles  of  three  characteristic  sizes  of  a  plasmoid  and 
three  velocity-  and  temperature-  components  in  the  plasmoid  center; 

-  time  profile  of  plasma  concentration  in  a  given  flow  point; 

-  coordinate  profiles  of  plasma  concentration  for  a  given  moment; 

-  line  of  equal  concentration  in  coordinate  planes  for  given  mo¬ 
ment  of  time. 

The  input  data  are  the  total  number  of  particles  in  a  plasmoid, 
initial  (t=0)  values  of  characteristic  sizes  and  their  spatial  deriva¬ 
tives,  and  also  temperature  given  in  the  origin  of  coordinates. 

Both  units  also  accumulate  data  arrays  in  separate  files  to  be 
further  used  for  i.g.  RF  refraction  /  scattering  calculations. 


2.2.  Representative  samples  of  calculations. 

Plasma  expansion  in  space  is  determined  by  either  the 
collision  processes  and  initial  conditions  or  inner  self  consis¬ 
tent  electric  fields.  These  fields  arise  due  to  high  electron  mobi¬ 
lity  and  are  proportional  to  local  density  and  temperature  gradi¬ 
ents.  In  rarefied  plasma  the  distribution  of  the  inner  fields  is 
strongly  influenced  by  the  thermal  conductivity  of  electrons  provi¬ 
ding  effective  transfer  of  energy  to  the  periphery  of  plasma  jet  or 
cloud  particularly  from  points  of  local  energy  releases.  The  ions 
are  accelerated  by  inner  plasma  fields  predominantly  in  the  direc¬ 
tion  of  the  maximal  gradient  (across  of  a  narrow  jet  or  along  the 
direction  of  heat  flux).  As  a  result  the  densities  of  energy  and 
substance  fluxes  in  expanding  plasmas  are  sufficiently  greater  than 
in  gases. 

However  these  peripheral  regions  are  extremely  important  sin¬ 
ce  they  determine  the  actual  boundaries  of  arising  inhomogeneities 
where  the  most  practically  interesting  processes  take  place  in 
fact.  For  instance  namely  these  regions  are  essentially  responsible 
for  RF  wave  refraction/scattering  and  for  electron  collection  by 
tether  plasma  contactors. 

Fig. 2. 1  displays  the  re¬ 
sults  of  calculation  for  free 
expansion  of  plasma  plume  of 
TSNIIMASH  made  experimental 
arcjet  plasma  source  EPICURE. 
The  procedure  of  input  data 
determination  based  on  labora¬ 
tory  measurements  and  more  de¬ 
tailed  comparison  of  the  lat- 
ters  with  theoretical  issues 


Fig.  2. 1 


obtained, is  presented  in  Sec. 3  It  is  believed  that  the  observed 
increase  of  curve  a(x)  inclination  at  low  values  of  x  represents 
the  adjacent  to  the  source  exit  region  of  one  -  two  exit  orifice 

diameters  long  where  energy 
exchange  between  plasma  speci¬ 
es  still  takes  place  before 
plasma  becomes  "frozen"  and 
stream  lines  become  straight. 

Fig. 2. 2  illustrates  jet 
expansion  of  very  high  Mach 
number  (M=14)  for  highly  focu¬ 
sed  charge  compensated  plasma 
Fig.  2.2  beam. 

The  characteristic  property  of  plasmoid  three-dimensional  ex¬ 
pansion,  namely  the  non-uniform  increase  of  its  sizes  clearly  seen 
in  Fig. 5  of  Appendix  A.  If  initial  velocities  are  low  compared  to 
the  mean  random  velocity,  vT~  [/T0/m,  the  resulting  velocity  is  ma¬ 
ximal  along  the  coordinate  axis  where  the  initial  size  of  a  plasmo¬ 
id  is  minimal.  Consequently  resulting  size  of  a  plasmoid  along  this 
coordinate  finally  occurs  to  be  maximal.  So  inversion  of  sizes  ta¬ 
kes  place. 

For  example,  a  plasmoid  initially  generated  as  a  "knitting 
needle"  stretched  along  axis  x,  later  transforms  in  a  disk  with  mi¬ 
nimal  x  dimension.  In  contrary  a  dense  disk  may  transform  in  a  long 
"knitting  needle".  The  reason  of  this  inversion  is  clear  enough: 
the  thermal  energy  initially  accumulated  by  a  plasmoid,  gradually 
transforms  in  kinetic  energy  of  expansion  predominantly  in  the  di¬ 
rection  of  maximal  gradient  of  electric  potential  or  in  other  words 
along  the  maximal  gradient  of  total  inner  pressure  V(nTi+nTe)  ~ 
VnTi+enVfl)  . 

The  results  of  calculations  were  further  used  for  comparison 
with  available  data  on  Plasma  free  expansion. 


3.  Comparison  of  results  obtained  with 
available  data  on  Plasma  and  Gas  free  Expansion 

3.1.  Comparison  with  Experimental  Data  In-House 

Thorough  measurements  of  plasma  free  expansion  characteristics 
were  undertaken  in  the  plume  of  the  experimental  arcjet  plasma  so¬ 
urce  EPICURE  developed  to  conduct  space  experiments  on  interaction 
between  RF-signal  and  an  electric  thruster  plume. 

The  device  basic  characteristics  are: 

-  input  power  _ <2  kW  (28V  x  70A); 

-  propellant  .  caesium; 

-  total  propellant  load _  0.3  kg; 

-  mass  flow  rate .  up  to  30- 10-6  kg/s. 

Two  operating  modes  were  chosen  for  comparison  of  experimental 
and  theoretical  data.  The  first  one  with  discharge  current  of  35  A 
occurred  to  be  corresponding  to  adiabatical  flow  with  criterion  n  < 
1  and  the  second  one  of  65  A  -  corresponding  to  isothermal  flow 
with  criterion  n  >  1. 

The  plasma  plume  basic  characteristics  -  electron  density,  n  , 
temperature,  T  ,  and  Mach  number,  M,  -  were  measured  by  Langmuir 
oriented  cylindrical  probes  located  along  the  vacuum  chamber  axis 
at  the  distances,  R,  of  0.3;  0.7;  and  1.2  m  from  the  plasma  source 
nozzle  exit  orifice.  Table  1  accumulates  the  plasma  parameters  mea¬ 
sured  at  the  flow  axis,  for  two  operating  modes  mentioned  above. 


Table  1 


operating  mode 

35  A 

65  A 

distance,  R,  m 

0.3 

0.7 

1.2 

0.3 

0.7 

1.2 

density,  n,  cm-3 

6  •  101 1 

1011 

2.5-1010 

2-1011 

3  •  101 0 

0.7-1010 

temperature,  T,  eV 

0.3 

- 

0.26 

1.8 

- 

1.6 

Mach  number,  M, 

2.5 

- 

2.6 

2.5 

- 

2.6 

The  plasma  source  was  mounted  on  a  rotating  gear  to  scan  the 
plasma  plume.  Fig. 3.2  shows  the  measured  angular  dimensionless  pro¬ 
files  of  n(0)/nc,  where  nc  is  the  value  of  electron  density  at  the 
flow  axis.  The  profile  at  R3  =  1.2  m  where  processes  of  energy  exc¬ 
hange  in  the  plume  are  believed  to  be  over  and  the  flow  becomes  co¬ 
nical,  was  used  to  determine  the  value  of  flow  divergency  factor, 
k=(a'0)"2,  and  particles  flow  rate,  N,  via  the  self-similar  angular 
distribution  for  the  operating  mode  of  35  A: 

ftvk  u3  1 

n(0)  - - — 

7TR3  2  u3  u  (l+k-tg20)p 


p 


1  + 


'XM2  /  ( M2  —  1 ) ; 


2 

U-DM32 

V=p-l/2 


[l-(a3/a)q] 

_  U-1)M2 
q  M2  -  1 


and  expression 

nc  =  n3(R3,0)  =  frvk/JtRc, 2  u3 . 

The  numerical  values  of  k  and  N  are  equal  3.0  and  3.4-1019  1/c 
correspondingly  so  that  the  criterion  n  is  equal  0.2,  e.i.  the  flow 
is  an  adiabatical  one  (II  <  1). 

These  numerical  values  of  k  and  N  were  used  as  the  input  data 
for  PC  program  DYNAMIC  to  calculate  the  profile  of  plasma  jet  cha¬ 
racteristic  transverse  size,  a(x)  along  the  jet  axis.  This  calcula¬ 
ted  profile  a(x)  is  plotted  in  Fig. 3.1  with  added  three  experimen¬ 
tal  points  revealing  good  agreement. 

The  dependence  a(x)  was  used  to  calculate  self-similar  angular 
distributions  according  to  relationships 

n  ac2  uc  (  R2  sin2 0^ "p 
nc  a2  u  V  a2  ) 

U-  .  t  +  -  fl  -  W  1 

uc2  M2  (y-1)  L  Va2  J  J 


Here  n=n(R,  0) ,  nc=n(R,0),  ac=a(R,  0),  a=a(R-cos0). 

The  calculated  electron  density  distributions  are  compared 


with  experimental  ones  in  Fig. 3. 2  for  three  distances  R=0.3,  0.7 

and  1.2m  from  the  nozzle  exit  orifice. 

Fig. 3.2  suggests  that  there  is  a  good  agreement  between  theory 
and  experiment  for  R=0.7  and  1.2m.  Comparison  of  theoretical  and 
experimental  distribution  at  R=0.3m  is  not  correct  since  the  ioni¬ 
zation  processes  are  not  still  over  within  this  nearest  downstream 
region.  This  suggestion  is  confirmed  by  spectroscopic  measurements. 

The  Robert's  approximation  (see  Ref.  [2]  from  report  [1])  are 
also  plotted  in  Fig. 3.2  revealing  agreement  with  the  experimental 
data  only  within  the  core  of  the  flow. 

For  the  operating  mode  of  65A  the  numerical  values  of  k  and  N 
are  found  to  be  5  and  4.2-1019l/c  correspondingly  so  that  the 
criterion  n  is  of  about  4,  e.i.  the  flow  is  actually  an  isothermal 
one  (II  >  1) . 

Using  the  same  procedure  to  determine  the  dependence  a(x)  (see 

Fig.  3.3)  again  the  self-similar  angular  density  distribution  was 

calculated  according  to  relationships 

n  ac2  uc  R2sin20^~p 

nc  a2  u  V  a2  ) 

u2  4  (&  \ 

—  =  1  +  -  In  | — | 

uc2  KM/-1  \3.„) 

The  resulting  curves  are  compared  with  experimental  electron 
density  distributions  in  Fig. 3.4  for  three  distances  R=o. 3,  0.7  and 
1.2m  from  the  nozzle  exit  orifice  also  revealing  better  agreement 
in  the  far  region  of  the  flow  field. 

Comparison  of  the  theory  and  experiment  presented  above,  con¬ 
firms  that  the  mathematical  model  of  plasma  free  expansion  correct¬ 
ly  describes  not  only  conical  flow  with  ultimate  velocity  and 
straight  flow  lines  in  far  region  but  also  in  adjacent  to  nozzle 
flow  region  where  plasma  heat  energy  accelerate  the  flow  and  flow 
lines  diverge. 
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3.2  Transversal  density  distribution 
in  adiabaticaly  expanding  plume 
of  Epicure  arcjet  (  n  <  1  ) 


o  experiment  (M^  2,5;  aj*  0,35) 

a(x) 


3.3  Isothermal  plume  expansion 
for  Epicure  arc  jet  (fl  >  1) 
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3.4  Transversal  density  distribution 
for  the  isothermal  mode  of  expansion 
of  Epicure  arcjet  plume  (  n  >  1  ) 


3.2.  Comparison  with  published  results 
3.2.1.  Traditional  thrusters 

Ref.  3  concerns  distributions  of  gas  density,  p,  and  substance 
and  impulse  fluxes,  pV  and  pV2,  in  Simulate  RSC  plumes.  Measured 
and  calculated  data  are  presented  as  angular  dependencies  in  polar 
coordinates  with  the  origin  at  the  center  of  the  nozzle  exit 
orifice.  The  measurements  were  made  at  relatively  large  constant 
distance  R  from  the  nozzle. 

Comparison  of  the  measurements  with  calculated  data,  based  on 
certain  mathematical  models  of  gas  flow,  suggests  that  experimental 
angular  distributions  of  gas  density  and  substance  and  impulse 
fluxes  are  more  sharp  in  the  vicinity  of  flow  axis  and  more  soft  at 
the  periphery  of  that  than  calculated  ones. 

In  the  self-similar  solution  for  adiabatically  expanding  jet 
the  transverse  distributions  of  the  velocity  head,  pV2,  does  not 
depend  on  angular  distribution  of  T (0)  and  V(0).  In  polar  coordina¬ 
tes  system  the  velocity  head  distribution  is  written  as 

mftumvk  uc  (l+tg20)2 

pV2  - - ; 

jtR2  um  (l+ktg20)1  +  p 

u  2 

—  =  1  +  b-  [  1  -  (1  +  tg20)q  ]  ; 

Urn  2 

Pm2  Mm2U-l)  2 

p  =  -  ;  q  =  -  ;  b  =  -  ; 

Mm2-1  Mm2-1  Mm2U-l) 

Here  um  and  Mm  are  maximal  velocity  and  Mach  number  at  the  flow 
axis,  0=0. 

It  is  taken  into  account  that 

a2  *  x2/k  =  R2cos20/k 

The  self-similar  distributions  plotted  in  Figs.  3  and  5  from 
Ref. 3  are  given  by 

pV2  (i+tg20)2  •  (i+b[l-(l+tg20)q]  }1/2 
pcVc2  (1  +  ktg20)1+p 


Using  function  f4  given  by 

rp  V  2  Un1/(1+p) 

fj  (0)  =  I*-2-1  (l+tg20)2  •— — |  =  1  +  ktg20 

LpV2  UmJ 

the  measured  angular  dependence  of  pV2  can  be  replotted  so  that  the 
inclination  of  resultant  straight  line  gives  the  value  of  flow 
divergency  factor,  k. 

Fig. 3.5  displays  plots  of  function  f±  (0)  for  N2  (^=1.4)  and 
mixture  N204/A-50  (^=1.31)  providing  values  k=7.9  and  k=7. 1  corres¬ 
pondingly.  For  these  two  working  substances  Fig  3.6  and  3.7  show 
derived  from  Ref.  3  theoretical  dependencies  pV27pcVc2  using  the 
same  coordinates. 

In  Fig  3.8  the  self-similar  angular  distributions  are  added  to 
those  derived  from  Ref [3]  with  corresponding  values  of  Mach  number. 
The  curves  of  pV  and  p  are  obtained  from  self-similar  angular 
distribution  (pV2)  using  experimental  dependence  V(0)  of  Fig  2  from 
Ref.  [3]. 

It  is  seen  that  the  self-similar  angular  dependencies,  displa¬ 
yed  in  Fig  3.6,  3.7  and  3.8  more  accurately  follow  the  experimental 
distribution  of  flow  parameters  than  simple  Robert's  approximation 
or  even  more  tedious  numerical  methods  generally  used.  Particular¬ 
ly,  the  self-similar  angular  dependencies  of  all  fluxes  are  more 
sharp  in  the  vicinity  of  flow  axis,  0  <  45°,  and  more  soft  at  its 
periphery,  0  >  60°,  than  those  calculated  using  other  flow  models 
of  freely  expanding  gas.  This  is  general  difference  between  experi¬ 
mental  and  numerical  results  discribed  in  [3]. 

Really,  the  expression  (pV2)  suggests  that  in  the  vicinity  of 
flow  axis,  20° <0<45° ,  e. i.  1/k  <  tg20  <1,  the  order  of  magnitude 
for  pV2  is  given  by  pV2  -  (ktg20)“ (p  +  1 }  whereas  at  flow  periphery, 
where  0  >  60°  and  tg20  »  1,  the  angular  dependence  of  pV2  is  more 
soft:  pV2  -  k" (p  +  1)  •  tg20"{p_1) . 
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.1  data  for  N2  and  N204/A-50  plum 
;]  plotted  In  special  coordinates 
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3.6  Experimental  and  theoretical  data  for  N2  plume 
from  Ref. [3]  plotted  In  special  coordinates 
In  comparison  with  self-similar  flow  model 
R  -  Roberts’  approximation 
S  -  Simons’  model 
A  -  Self-similar  model 


M  =3,95 
K  =7.9 


3.7  Experimental  and  theoretical  data  for  N204/A-50 
plume  from  Ref.  [3]  plotted  in  special  coordinates 
in  comparison  with  self-similar  flow  model 
R  -  Roberts’  approximation 
S  -  Simons’  model 
A  -  Self-similar  model 
Inv  -  Inviscid 

VT.VL  -  Viscid  turbulent  and  laminar 
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ANGLE  FROM  CENTERLINE,  DEG 

Fig.  2  Normalized  mass  flux  and  velocity  measurements  in  the  plume 
of  a  1/10-scale  Lunar  Module  RCS  engine. 


SIMONS'  MODEL  (<UBL>/Umas  *  0.68) 


Fig.  4  Comparison  of  measured  density  ratio  with  predictions  from 
Simons’  model. 


Fig.  3  Comparison  of  measured  dynamic  pressure  profile  with  Boyn¬ 
ton’s  predictions. 


Fig.  5  Comparison  of  Chirivella’s  measured  flux  ratios  with  predic¬ 
ted  density  ratio  from  Simons’  model. 


3.8  Experimental  and  theoretical  plots  of  Ref. [3J 
with  self-similar  curves  added  for  comparison 


3.2.2.  Electric  Thrusters  Plumes 


Operation  of  onboard  Stationary  plasma  (SPT)  or  Anode  layer 
(TAL)  electric  thrusters  results  in  arising  in  space  of  expanding 
ion  beam  with  the  following  representative  initial  parameters  : 
ion  energy  E0=  rnvm2/2  =  270  eV; 

electron  temperature  T  ~  1  *  5  eV; 

values  of  square  Mach  number  M2  ~  mv02/T0Y  ~  300  +  60; 
corresponding  to  criterion  n  =  6Te7taa Ve2N  ~  20  +  1000. 
Consequently  ion  beam  expansion  is  isothermal  under  the 
condition  of  M2  -»  °°. 

Taking  into  account  that  for  the  nearest  flow  region 
a0  <2  =  const  =  1/k,  for  a/a0  «  exp(mv02/2T0 ) 
we  obtain  conical  beam  expansion. 

Ref.  4  presents  the  distribution  of  ion  beam  current  density 
j(x,  r)  measured  in  cylindrical  coordinates.  There  is  a  self-similar 
relationship  for  this  quantity: 

eNk  1 

2t tx2  (l+kr2/x2  )3/2 

Theoretical  and  experimental  dependencies  can  be  compared  using 
two  independent  methods.  The  ion  beam  current  density  at  the  axis 
is  expressed  via  total  ion  beam  current,  1^  =  eN,  and  factor  of  flow 
divergency,  k. 

eNk 


The  measured  transverse  dependence  of  jx  can  be  replotted  in 
the  special  coordinates 

f2(if\2)  =  ( Jc /Jx ) 2 /3  =  1  +  k-r2/x2 

so  that  the  inclination  of  resultant  straight  line  gives  the  value 
of  flow  divergency  factor,  k. 

Function  f2  is  plotted  in  Fig.  3.9  for  a  few  transverse  secti- 


ons  of  the  ion  beam  for  x  ranging  from  0.52  to  3.7  m.  Fig. 3.9  disp¬ 
lays  jc (x)  curve  calculated  for  k  =9.5  and  experimental  points  for 
axial  values  of  jc.  Experimental  and  calculated  lines  of  equal  va¬ 
lues  of  J'x  are  plotted  in  Fig.  3.9. 

All  three  Figures  suggest  that  experimental  and  calculated  da¬ 
ta  satisfactorily  correlate  for  x  <  3  m.  Under  the  conditions  on 
ambient  pressure  for  experiments  of  Ref.  4,  the  charge  exchange 
free  path  length  for  Xe+  ions  is  evaluated  on  the  level  of  about 
3m.  Therefore  the  observed  difference  between  experimental  and  cal¬ 
culated  data  for  x  >  3  m,  is  believed  to  be  due  to  ion  beam  scatte¬ 
ring  in  the  ambient  gas. 

In  Ref. [5]  the  angular  distributions  of  ion  beam  current  densi¬ 
ty,  j (0) ,  and  ion  energy,  E(0),  for  SPT-70  electric  thruster  plume 
are  measured  in  polar  coordinates  at  constant  distance  R,  from  the 
exit  orifice  plane.  The  self-similar  solution  gives  the  following 
angular  dependencies  for  these  quantities: 

efe-[/E/2m  (l+tg20)2 
2jtR2  (l+ktg20)2’ 

Again  the  appropriate  function  f3  can  be  plotted  in  the  speci¬ 
al  coordinates  so  that  the  inclination  of  resultant  straight  line 
gives  the  value  of  flow  divergency  factor,  k.  From  Fig. 3.10  the  re¬ 
sultant  value  of  k  =  9.6,  practically  is  equal  to  that,  k  =  9.5, 
determined  for  another  SPT  version  of  the  same  design  [4].  Although 
in  Refs.  4  and  5  the  SPT  electric  thrusters  of  the  same  design  were 
studied,  the  measuring  methods  were  principally  different  and  vacu¬ 
um  facility  conditions  also  differed.  Coincidence  of  laboratory  me¬ 
asurements  processing  results  for  these  different  cases  suggests 
that  the  used  analytical  method  is  correct  enough. 

Let  us  compare  the  experimental  data  on  ion  beam  expansion  repor¬ 
ted  in  Refs. 6, 7  and  8  for  SPT-100  and  TAL  D-55  electric  thrusters 
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3.9  Comparison  of  experimental  data  from  Ref. [4] 
for  the  SPT  type  electric  thruster  exhaust 
beam  with  the  self-similar  flow  model  of 
isothermal  expansion  mode 

-  self-similar  model 

-  Robert's  approximation 


3.10  Comparison  of  experimental  data  from  Ref. [5-8] 
for  SPT-70,  SPT-100  and  TAL  D-55  electric  engine 
plumes  with  the  self-similar  flow  model  of 
Isothermal  expansion  mode 
•  SPT-70  experiment  [5] 

0  SPT-100  experiment  [6,7] 

x  TAL  D-55  experiment  [6,8] 

-  self-similar  model 


of  different  design  under  the  conditions  of  the  same  operating  mo¬ 
de:  accelerating  voltage  of  300  V,  discharge  current  of  4.5  A  and 
equal  Xe  mass  flow  rate.  Value  of  electron  temperature,  Te,  is  ar¬ 
bitrary  taken  at  its  upper  limit  of  5  eV. 

Self-similar  ion  beam  current  density  is  given  by; 

efr-k  1 

Ji  =  env  = - : 

2jiR2cos30  (l+ktg20)3/S 

Again  the  appropriate  function  f4  being  plotted  in  the  special 
coordinates,  gives  the  value  of  flow  divergency  factor,  k.  Fig. 3. 10 
shows  that  the  resultant  values  of  k  occur  to  be  equal  26  both. 
This  suggests  that  two  ion  accelerators  using  almost  the  same  phy¬ 
sical  principles  despite  of  different  design,  provides  practically 
the  same  beam  focusing  with  characteristic  divergence 

arctg(a')  =  11° 

It  is  also  shown  that  plumes  of  Stationary  plasma  (SPT)  or 
Anode  layer  (TAL)  electric  thrusters  producing  high  velocity  ion 
beam,  are  characterized  with  more  spread  substance  fluxes  than  gas 
flows  due  to  action  of  radial  electric  fields  produced  by  the 
electronic  specie. 


Discussion 


The  demonstrated  above  good  agreement  between  experimental  and 
theoretical  data  suggests  that  here  presented  approach  to  formula¬ 
tion  and  solution  of  plasma  free  expansion  dynamical  problem  is 
fruitful  and  correct  enough. 

This  approach  can  be  further  developed  to  evaluate  interaction 
of  freely  expanding  plasma  with  environment  both  in  vacuum  facili¬ 
ties  and  space. 

In  Ref.  [7]  was  noted  that  lowering  of  background  pressure  in 
vacuum  chamber  results  in  electric  thruster  plume  increase  in  cont¬ 
radiction  to  predictions  based  on  mechanisms  of  beam  scattering  in 
residual  atmosphere.  It  means  that  electric  thruster  exhaust  beam 
divergence  in  deep  vacuum  of  space  flight  may  be  significantly  gre¬ 
ater  than  that  on-ground  measured. 

Authors  of  Ref. [93  attempted  to  explain  temporal  degradation 
of  SPT  operating  characteristics,  taking  into  account  processes  of 
plume  interaction  with  ambient  gas  and  subsequent  impact  to  accele¬ 
rating  mechanisms  in  SPT. 

The  described  above  method  of  plasma  jet  parameters  calculati¬ 
on  also  may  be  used  to  study  the  mechanisms  of  plasma  jet  interac¬ 
tion  with  ambient  gas  resulting  in  correction  of  on-ground  measured 
parameters  of  electric  thrusters  for  the  actual  conditions  of  space 
environment.  Particularly,  the  unexpected  results  of  Ref. [7]  may  be 
satisfactorily  explained. 

This  approach  also  allow  to  evaluate  the  expected  impact  of 
the  spacecraft  ambient  atmosphere  including  return  fluxes  of  subs¬ 
tance  released. 

The  formulation  and  solution  of  dynamical  problem  on  plasma 
free  expansion  presented  in  Report  [ i] ,  may  be  readily  updated  to 


take  into  account  electric  currents  flowing  in  plasma  and  ambient 
media,  e.g.  in  metallic  walls  of  vacuum  chambers  and  in  ionosphere. 
The  former  updating  could  be  further  used  to  improve  correctness  of 
on-ground  measurements  and  the  latter  one  is  perspective  for  deve¬ 
lopment  of  theory  of  plasma  contactors  for  electrodynamical  tethers 
in  space. 

The  solutions  of  plasma  expansion  problem  obtained,  are  valid 
in  relatively  near  region  of  plasma  flow  where  forces  of  interacti¬ 
on  of  flowing  plasma  and  magnetic  field  are  low  compared  to  pressu¬ 
re  gradients. 

The  relation  between  these  forces  is  determined  by  so  called 
parameter  of  magnetic  interaction, 

6u2  B2  jra3  a  ‘ 

S  =  — - 

NT 

The  characteristic  size  of  plasma  flow,  where  S  <  1  and  action 

of  magnetic  field  is  relatively  insufficient  is  given  by 

/'NTk2  a  1/3 

XB  =  - 

v6u2b2; 

For  example  this  estimated  size  for  plume  of  Epicure  arcjet 
(see  Section  3.1),  operating  in  space  is  of  about  3, 2  m  with  due 
regard  for  action  of  the  Earth  magnetic  field.  For  plume  of  elect¬ 
ric  thruster  of  SPT  type  (see  section  3.2.2)  this  size  is  only  2  m. 

Therefore  for  many  real  situations  the  action  of  the  Earth 
magnetic  field  must  be  taken  into  account  while  evaluating  interac¬ 
tion  of  on-board  devices  plasma  plumes  with  environment  and  spacec¬ 
raft  structure,  and  calculating  RF  scattering  effects. 


Conclusion 


The  dynamical  problem  describing  evolution  of  both  fully  ioni¬ 
zed  stationary  supersonic  plasma  jet  and  plasmoid  injected  from  a 
spacecraft  is  formulated  and  solved  in  a  self-similar  form. 

Because  of  the  self-similar  approach  and  reasonable  simplifi¬ 
cation  the  initial  combined  partial  differential  equations  of  two 
components  plasma  dynamics  were  rearranged  in  a  set  of  separable 
equations  so  that  the  problem  was  reduced  to  quadratures  and  to  so¬ 
lution  of  ordinary  nonlinear  differential  equations. 

Analytic  forms  of  the  solution  were  obtained  for  practically 
interesting  ultimate  thermal  modes  of  plasma  flows  including 

-  hypersonic  adiabatic  plasma  or  gas  jet  expansion; 

-  isothermal  expansion  of  highly  accelerated  ion  beam  and 

-  ellipsoidal  plasmoid. 

The  dynamical  problem  solution  obtained,  properly  describes 
either  the  adjacent  to  a  nozzle  region  of  non-linearly  expanding 
flow  or  the  far  region  of  inertial  expansion. 

The  most  noticeable  feature  of  three  dimensional  plasmoid  ex¬ 
pansion  is  the  inversion  of  its  shape. 

Validity  of  theoretical  result  obtained  is  confirmed  by  the 
way  of  comparison  with  in-house  experimental  and/or  published  data. 
Good  agreement  between  experimental  and  theoretical  data  is  demons¬ 
trated  for  traditional  and  electric  thrusters  plumes  both.  It  is 
also  shown  that  presented  mathematical  model  of  freely  expanding 
plasma  jet  more  precisely  follows  the  experimental  data  than  Ro¬ 
bert's  and  Narasimha's  approximations  often  used. 

The  obtained  dynamical  problems  solutions  were  put  in  the  base 
of  PC  Program  DYNAMIC  developed  to  calculate  characteristics  of 
free  expansion  of  artificial  plasma  formations  in  space.  Compact 
forms  of  solutions  used,  result  in  great  saving  of  run  time.  The 
Program  DYNAMIC  allows  "to  visualise"  plasma  flow  field  and  provi¬ 
des  data  for  numerous  application  including  i.g.  RF  refraction  / 
scattering  calculations,  necessary  for  determination  of  communica¬ 
tion,  navigation  and  ranging  conditions  for  a  spacecraft. 


*  *  * 
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Appendix  A 


Program  for  calculation 
of  plasma  free  expansion  in  vacuum 
DYNAMIC 


Program  Description 
Content: 

1.  Assignment,  content  and  conditions  of  use 

2.  How  to  run  the  Program 

3.  Unit  JET 

3.1.  Physical  base 

3.2  General  Algorithm 

3.3  Input  data 

3.4  Input  data  ranges 

3.5  Results  description 

4.  Unit  PLASMOID. 

4.1.  Physical  base 

4.2  General  Algorithm 

4. 3  Input  data 

4. 4  Input  data  ranges 

4.5  Results  description 


1.  Assignment,  content  and  conditions  of  use 


1.1.  Program  DYNAMIC  is  intended  to  calculate  characteristics  of 
free  expansion  of  artificial  plasma  formations  (APF)  arising  in  space 
due  to  operation  of  onboard  plasma  sources. 

1.2.  It  consists  of  two  units  -  JET  and  PLASMOID. 

The  former  calculates  characteristics  of  plasma  Jet  free  expansi¬ 
on  in  vacuum  for  plasma  sources  of  steady  state  mode  operation. 

The  latter  calculates  characteristics  of  plasmoid  free  expansion 
in  vacuum  for  plasma  sources  of  pulse  mode  operation. 

1.3.  Program  DYNAMIC  provides  the  following  as  a  result  of  calcu¬ 
lation: 

-APF  dimensional  parameters; 

-spatial  and  temporal  distributions  of  plasma  concentration  and 
temperature; 

-configuration  of  equal  plasma  concentration  surfaces  . 

The  Program  presents  the  results  of  calculation  run  graphically 
on  a  display  and  also  accumulates  resultant  data  arrays  in  separate 
files  of  current  directory  for  further  use. 

1.4.  To  put  DYNAMIC  run  You  need  a  IBM  compatible  computer  with  a 
keyboard  of  101  keys  and  more  than  100  kilobytes  memory  available  in 
current  directory. 

1.5.  Program  DYNAMIC  is  developed  in  Turbo  Pascal  7.0  integrated 
environment  to  be  run  with  Microsoft  MS  DOS  6.2  or  later  version  or 
compatible  one. 

Program  DYNAMIC  contains  the  following  four  files: 

-DYNAMIC. EXE  (the  main  file); 

-STRM-IN  (input  data  file  of  test  example  for  JET  unit); 

-PLSD-IN  (input  data  file  of  test  example  for  PLASMOID  unit); 

-READ. ME  (text  file  of  user  manual). 

1.6.  Program  DYNAMIC  provides  the  following: 

-choice  of  calculation  version; 

-data  entry,  input  data  correction; 

-calculation  runs  and  resultant  data  soft  copy; 

1.7  The  run  time  of  one  calculation  version  is  of  about  1  minute 
for  RC  AT  486. 


2.  How  to  run  the  Program 


2.1.  Getting  Program  DYNAMIC  run  Is  simple.  Insert  DYNAMIC  floppy 
disk  delivered  In  an  appropriate  drive  or  Enter  the  directory  on  hard 
disk  where  You  have  stored  DYNAMIC  files  to  save  original  floppy  disk 
and  choose  DYNAMIC. EXE.  resulting  in  appearance  of  window  1  (see 
Fig. 1)  with  brief  description  of  Program  DYNAMIC. 
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Fig.  1 

2.2.  Press  ENTER  to  display  Program  DYNAMIC  MENU  (see  Fig. 2).  You 
can  use  the  arrow  keys  to  choose  one  of  two  Program  units  or  press  ESC 
to  interrupt  the  Program  run. 

2.3.  If  You  made  Your  choice,  press  ENTER  to  open  window  3  (see 
Fig.  3  ).  In  this  window  You  can  type  new  input  data  for  calculation 
pressing  ENTER  after  each  typing  or  choose  ESC  to  start  execution  of 
the  test  example.  For  more  information  about  input  data  entry,  see 
descriptions  of  Program  units. 


Fig. 2 


Initial  dimensionless  parameters 


WARNING!  Window  3  displays  input  data  for  test  example  contained 
in  *-in  file  of  the  floppy  disk  delivered.  While  entering  new  input 
data  in  window  3  You  overwrite  old  ones.  Column  "New  Data"  must  be  ty¬ 
ped  down  all.  By  choosing  ESC  You  get  run  test  version  with  updated 
input  data  which  do  not  correspond  to  information  of  column  "Test"  in 
the  table  of  window  3.  Therefore  You  may  save  time  if  You  edit  *-lN 
file  before  getting  DYNAMIC. EXE  run.  In  this  case  while  being  in  win¬ 
dow  3  You  merely  press  ESC  to  perform  test  version  with  already  upda¬ 
ted  input  data. 

WARNING!  If  current  directory  does  not  contain  *-lN  file,  press 
ENTER  and  write  input  data  from  the  screen. 

2.4.  The  result  of  Program  DYNAMIC  run  appears  in  window  4.  It 
looks  like  Fig. 4  for  unit  JET  and  like  Fig. 5  for  unit  PLASMOID.  For 
more  information  about  representation  of  the  Program  run  results,  see 
descriptions  of  Program  units. 

2.5.  As  the  Program  run  is  over,  press  ENTER  to  return  to  MENU  in 
window  2.  Here  You  may  choose  the  new  calculation  version  or  escape 
Program  DYNAMIC  by  pressing  ESC. 


3.  Unit  JET 


3.1.  Physical  base 

JET  unit  calculates  characteristics  of  plasma  jet  free  expansion 
in  vacuum  for  plasma  sources  of  steady  state  mode  operation. 

The  general  assumption  is  that  two  components  plasma  is  fully  Io¬ 
nized  and  Its  flow  Is  adiabatic  one. 

Self-consistent  plasma  flow  in  a  steady  state  axial  symmetric  jet 
Is  described  with  the  set  of  equations  using  cylindrical  coordinates  x 
and  r: 

dnu  +  drnu  =  0; 

3x  rdr 

(  du  du  \  dnT 

mn  |u —  +  v —  |  = - ; 

V  dx  dr  )  dx 

f  dv  dv  \  dnT 

mn  |u —  +  v —  I  = - ; 

V  dx  dr  J  dr 

n  (  dT  dT  \  (du.  drv  a 

-  |u —  +  v —  |  +  nT| —  +  -  I  =  Vq  , 

Y~1  V  dx  dr  J  Vdx  dr  J 

The  self-similar  form  of  solution  for  the  system  (4)  Is  written 

as: 

u  =  uc  (x)  •  y (ti) ;  v  =  u-da/dx-ti;  T  =  tc(x)*t(ti);  nu  =  Nv  /Jta2-f(t\), 
r  =  a(x)  is  a  certain  flow  current  line  (see  Fig.  4). 

For  three  functions  uc(x),  Tc (x)  and  a(x)  there  are  four  equations: 

muc  2  'X  muc  a 

—  +  —  Tc  =  const;  -  (uca’)'  =  Ct; 

2  *-l  Tc 

T 

Tc  i  /  (^_ i )  •  uc a2  =  const;  —  a{^_1)  =  const; 

uc 

where  C1  and  C2  are  separation  constants  for  variables  x  and  t\ 
correspondingly. 

The  latter  Is  Identical  to  the  adiabatic  equation  when  uc  const 
and  C2  =  2%.  For  adiabatic  expansion  the  flow  velocity  always  reaches 
its  maximal  value  and  remains  practically  constant,  i.e. 


u2  -*  um2  =  u02  [1  +  2 -Mo" 2  •  (K  -if"1  ] . 

With  C* =C2 ,  we  obtain  the  following  relationships  for  jet  expansi¬ 
on  along  the  axis: 

uc  =  um  ;  Tca2^_1)  =  T0a02(?Ll)  ;  rnumaa"  =  2T0  (a0/a)2 (f~1 )  (1) 

For  rarefied  plasma  flows  with  uniform  temperature  in  the  flow  cross 
sectional  plane  in  vicinity  of  the  flow  pole  due  to  high  thermal  con¬ 
ductivity  of  the  plasma  9T/9r=0,  Tc=Tc(x),  the  problem  solution  is 
written  as: 

_  N-  0-1/2)  1  _  _  um 

jra2um  (l+r2/a2)  *  (l+r2/a2)1/2 

T  =  Tc0  •  (a0/a)2  } ;  a  =  a0  +  am  (x-x0)  ;  (2) 

am2  =  a02  +  2  /  (Of-l)-M02)  ; 

These  formulas  suggest  that  the  distributions  of  plasma  flow  pa¬ 
rameters  in  the  far  region  of  the  jet  where  a  »  a0  are  defined  by  the 
particles  flow  rate,  N,  the  ultimate  plasma  velocity,  um,  and  the  fac¬ 
tor  of  flow  divergency,  k  =  l/(a’m)2. 


3.2  General  Algorithm 


The  first  step  Is  calculation  of  profile  of  plasma  jet  characte¬ 
ristic  transverse  size  ,a(x),  along  the  jet  axis,  x.  In  the  way  of 
Kutta-Merson  Integration  of  the  ordinary  nonlinear  differential  equa¬ 
tion  (1).  The  Input  data  are  given  at  x=0:  a0  =1,  T0=l.  Values  of  'K, 
a0 '  h  M2  are  specified  by  a  user.  A  user  also  chooses  mutually  corre¬ 
lated  Integration  step  h  and  area  s  taking  Into  account  available  PC 
capacity  and  display  resolution.  For  Instance  If  You  specify  s=10,  the 
proper  choice  for  h  Is  0.01  and  0.1  for  s=100  to  obtain  curves  capital 
enough  on  a  screen  and  to  save  time  run.  The  resultant  values  of  a(x) 
at  every  Integration  step  are  collected  In  file  A. OUT  and  correspon¬ 
ding  values  of  x  coordinate  are  stored  In  file  S.OUT.  These  data  are 
reserved  for  further  use  after  Program  run  Is  over. 

Integration  run  Is  stopped  as  condition  x>s  Is  satisfied  and  fi¬ 
les  A.  OUT  and  S.OUT.  become  open  for  readout  and  use  In  calculation  of 
spatial  distributions  of  plasma  concentration  n(x)  and  temperature 
T(x)  using  expressions  (2). 

The  expression  for  n(x)  of  set  (2)  Is  also  used  to  calculate  con¬ 
figuration  of  equal  plasma  concentration  lines  via  variation  of  trans¬ 
verse  coordinate  r. 


3.3  Input  data 

The  calculations  are  run  In  the  dimensionless  form.  The  scale  of 
linear  dimensions  Is  initial  transverse  size  of  the  jet,  a0,  and  the 
temperature  and  plasma  concentration  scales  are  Initial  values  T0  and 
n0. 

Input  parameters  for  JET  unit  are  the  following: 

g  =  specific  heat  ratio  'K; 


M  =  Mach  number; 

h  =  integration  step; 

s  =  integration  area  along  the  Jet  axis; 

n/n0  =  equal  plasma  concentration  line  parameter; 

a'  =  flow  divergence. 


3.4  Input  data  ranges 


Identifier 

Range 

Note 

8 

1.  .  .5/3 

=  5/3  for  test  example 

M 

2.  .  .  15 

=2.5  for  test  example 

h 

0.001...1 

choice  of  correlation 

s 

1. . . 100 

between  h  and  s  see  3.2 

n 

0.999.  .  .0.001 

=0.8  and  0.6  for  test  example 

a' 

0. . .  1 

=0.2  for  test  example 

3.5  Results  description 

Results  of  JET  unit  run  are  profile  of  plasma  Jet  characteristic 
transverse  size  ,  a  ,  along  the  Jet  axis,  x,  spatial  distributions  of 
plasma  concentration  n(x)  and  temperature  T(x)  and  equal  plasma  con¬ 
centration  lines  for  two  levels  (0.8  and  0.6  for  test  example).  Screen 
view  is  shown  in  Fig.  4.  The  additional  output  information  are  data 
saved  in  files  A. OUT  and  S.OUT. 


Fig. 4 


4. Unit  PLASMOID. 


4.1.  Physical  base 

Unit  PLASMOID  calculates  characteristics  of  plasmoid  free  expan¬ 
sion  in  vacuum  for  plasma  sources  of  pulse  mode  operation. 

For  three-dimensional  non-steady-state  plasmoid  the  system  of  dy¬ 
namic  equations  in  Cartesian  coordinates  is  written  as: 
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Here  T=Te+Ti  is  the  sum  of  ion  and  electron  temperatures.  Assumed 
dependencies  of  velocities  on  coordinates  are  linear: 


x 


u  =  a 


v  =  b 


y 


w 


.  z 

c  — 


a  b  c 

where  a  =  a(t),  b  =  b(t)  c  =  c(t)  are  characteristic  sizes  of  the 
plasmoid  along  axes  x,  y,  z  respectively.  Spatial  distribution  of  plasma 
concentration  and  temperature  within  a  plasmoid  is  written  as: 

N  v 

—  ■  h  (n)f2  (Of3  (t)  ; 

8abc  (3) 


n  = 


T  =  Tc  (x)  x  (  n,  t  x,  ) 

where  v  is  a  normalization  factor  defined  from  the  integral 

00 

Iff  n  dxdydz  =  N, 

00 

where  N  is  the  total  amount  of  particles  in  the  plasmoid,  and  f1# 

x  y  z 

f2  and  f3  are  functions  of  self-similar  variables  ri= —  ;£= —  ;t,= — 

a  b  c 


such  that  f±  (0)  =  f2  (0)  =  f3  (0)  =  1. 

The  representative  solution  corresponds  to  the  case  of  uniform 

temperature  throughout  the  plasmoid  e.i.  t  =  1  and  Vq  =  0.  This  case 

is  realized  for  "hot"  plasmoid  when  the  role  of  heat  conductivity  in 

plasma  is  sufficient,  especially  under  condition  of  no  energy  exchange 

between  the  plasmoid  and  the  environment.  In  this  case: 

r  x2  )  f  y2  )  (  z*  \ 

fi  =  exp - - C±  | ;  f2  =  exp - -C2  ;  f3  =  exp - -C3  ;  (4) 

V  2a2  J  V  2b2  J  V  2c2  J 

Characteristic  sizes  of  a,  b,  c  are  specified  as  those  ones  correspon¬ 
ding  to  "knee"  in  dependence  of  n  on  coordinates  and  to  the  condition: 

f/'  (1)  =  f2"  (1)  =  f3"  (1)  =  0 


The  separation  constants  can  be  also  determined  from  these 
condition: 

Ci  =  C2  =  C3  =  1 

For  a  "heat  insulated"  plasmoid  with  uniform  temperature  the 
energy  equation  (3)  is  reduced  to  the  adiabatic  relationship: 


m  f  n  f 1  ra0b0c0f 

Tc  =  T0  - 1  =  T0  - - 


V  n0J 


(5) 


label 

For  definition  of  time  dependent  sizes  of  the  plasmoid  there  is  a 
set  of  non-linear  differential  equations  : 

f-: 


maa  =  mbb  =  mcc  =  To 


pk,  b0  c0  j  1 
V  a  b  c  J 


(6) 


Here  dots  denote  derivatives  with  respect  to  time,  t.  The  charac¬ 
teristic  property  of  this  three-dimensional  expansion  is  non-uniform 
increase  of  plasmoid  sizes.  If  initial  velocities  are  low  compared  to 
the  mean  random  velocity,  vT~  |/f0/m,  the  resulting  velocity  is  maximal 
along  the  coordinate  where  the  initial  size  of  a  plasmoid  is  minimal. 
Consequently  resulting  size  of  a  plasmoid  along  this  coordinate  final¬ 
ly  occurs  to  be  maximal.  So  inversion  of  sizes  takes  place. 

For  example,  a  plasmoid  initially  generated  as  a  "knitting  need- 


le"  stretched  along  axis  x,  later  transforms  In  a  disk  with  minimal  x 
dimension.  In  contrary  a  dense  disk  may  transform  in  a  long  "knitting 
needle".  The  reason  of  this  inversion  is  clear  enough:  the  thermal 
energy  initially  accumulated  by  a  plasmoid,  gradually  transforms  in 
kinetic  energy  of  expansion  predominantly  in  the  direction  of  maximal 
gradient  of  electric  potential  or  in  other  words  along  the  maximal 
gradient  of  total  inner  pressure  VtnTi+nTg)  ~  VnTi+enVO  . 


4.2  General  Algorithm 


The  first  step  Is  calculation  of  temporal  dependencies  of  plasmoid 
linear  dimensions  a(t),  b(t)  and  c(t)  in  the  way  of  Kutta-Merson  in¬ 
tegration  of  combined  ordinary  nonlinear  differential  equations  (6). 
The  input  data  are  given  at  t=0:  T0=l,  m0=l,  a0 ‘=b0 ‘=c0 ‘=0.  Values  of 
Y,  M,  a0,  b0,  c0  are  specified  by  a  user.  A  user  also  chooses  mutually 
correlated  integration  step  h  and  temporal  area  t  taking  into  account 
available  PC  capacity  and  display  resolution.  For  instance  if  You  spe¬ 
cify  t=10,  the  proper  choice  for  h  is  0.01  and  0.1  for  t=100  to  obtain 
curves  capital  enough  on  a  screen  and  to  save  time  run.  The  resultant 
values  of  a(t),  b(t)  and  c(t)  as  well  as  derivatives  a'(t),  b'(t)  and 
c'(t)  at  every  integration  step  are  collected  in  files  A. OUT  B. OUT  and 
C. OUT,  A1.0UT  Bl. OUT  and  Cl. OUT  correspondingly  whereas  values  of  "t 
coordinate"  are  stored  in  file  T.OUT.  These  data  are  reserved  for 
further  use  after  Program  run  is  over. 

Integration  run  is  stopped  as  condition  t>x  is  satisfied  and  fi¬ 
les  A. OUT  B. OUT  and  C. OUT,  A1.0UT  B1.0UT  and  Cl. OUT  become  open  for 
readout  and  use  in  calculation  of  temporal  distributions  of  plasma 
temperature  T(t)  and  concentration  n(t)  at  the  plasmoid  center  using 
expressions  (5)  and  (3)  combined  with  (4)  correspondingly  . 

The  expression  (3)  combined  with  (4)  is  also  used  to  calculate 
configuration  of  equal  plasma  concentration  lines  of  n/n0  =0.5  level 
for  two  fixed  moments  t^SO-h  and  t2=500-h. 

4. 3  Input  data 

The  calculations  are  run  in  the  dimensionless  form.  The  scale  of 
linear  dimensions  is  initial  maximal  size  of  the  plasmoid,  a0.  Gene¬ 
rally  the  plasmoid  is  characterized  with  three  initial  sizes  -  maxi¬ 
mal,  mean  and  minimal.  In  the  test  example  the  plasmoid  is  initially 


an  elliptic  cylinder  with  the  maximal  size  being  tenwise  greater  than 
the  minimal  one.  The  temperature  and  plasma  concentration  scales  are 
initial  values,  T0  and  n0. Values  of  t0=i/ma02/T0  is  the  temporal  scale. 
Input  parameters  for  PLASMOID  unit  are  the  following: 
aO  =  maximal  initial  size; 
bO  =  mean  initial  size; 
cO  =  minimal  initial  size; 
g  =  specific  heat  ratio 
h  =  integration  step; 
t  =  duration  of  free  expansion. 


4.4  Input  data  ranges 


Identifier 

Range 

Note 

aO 

0. . .  1 

=  1.0  for  test  example 

bO 

0. . .  1 

=0.2  for  test  example 

cO 

0. . .  1 

=  0.02  for  test  example 

h 

O 

o 

o 

1-^ 

choice  of  correlation 

t 

1. . . 100 

between  h  and  t  see  4.2 

g 

1. .  .5/3 

=  5/3  for  test  example 

4.5  Results  description 

Results  of  PLASMOID  unit  run  are  the  following  (screen  view  is 
shown  in  Fig. 5)  : 

-  time  profiles  of  three  characteristic  sizes  and  three  velocity 
components  of  the  plasmoid;  ('Fig. 1') 

-  time  profile  of  plasma  concentration  and  temperature  in  the 
center  of  the  plasmoid; ('Fig. 2’) 

-  coordinate  profiles  of  plasma  concentration  for  moment  t2: 


-  line  of  equal  concentration  in  coordinate  planes  for  moments 
tj  and  t2  -  n/nc=0.5.  (’Fig.  4, 5’) 

The  additional  output  information  are  data  saved  in  files  *.0UT. 
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Appendix  B 


Proposals 
on  further  stages 

of  Development  of  dynamical  and  mathematical  models 
of  plasma  flows  in  space 

with  due  regard  for  action  the  Earth  magnetic  field. 

This  Appendix  submits  a  set  of  new  proposals  for  EOARD  that 
could  be  a  perspective  of  future  cooperation.  This  set  encloses  a 
part  of  our  professional  interests  where  we  have  a  sufficient  theo¬ 
retical  and  experimental  background. 

The  proposals  submitted  are  based  on  the  main  principle  that 
every  given  proposal  covers  one  or  a  number  solved  problem  as  a  fi¬ 
nal  result  available  for  a  user.  As  a  rule  validity  of  this  result 
is  evaluated  by  the  way  of  comparison  with  in-house  experimental 
and/or  published  data. 

The  set  of  given  new  proposals  concerns  three  types  of  artifi¬ 
cial  plasma  formations  (APFs)  which  can  arise  in  space  particularly 
in  vicinity  of  a  spacecraft  as  a  typical  result  of  normal  operation 
of  certain  onboard  systems  or  may  be  intently  created  for  scienti¬ 
fic  or  applied  goals:  stationary  plasma  plumes,  plasmoids  or  ioni¬ 
zed  clouds  and  gas-plasma  formations  as  a  result  of  natural  ioniza¬ 
tion  of  neutral  gas  releases.  Adequate  description  and  parameters 
calculation  of  every  APF  type  presents  an  independent  chain  of  dy¬ 
namical  problems,  methods  and  codes.  Every  chain  is  multi-link  sin¬ 
ce  for  every  APF  type  one  may  define  a  few  (2  -  4)  stages  of  expan¬ 
sion  in  dependence  on  predominant  mechanisms  and  degree  of  interac¬ 
tion  of  injected  substance  with  main  dynamical  factors  of  space 
flight,  namely  the  Earth  magnetic  field  and  environment. 


The  first  stage  corresponds  to  non-perturbed  expansion  of  dense 
plasma  flow  core  where  the  plasma  density  and  pressure  are  much 
greater  than  the  density  of  environmental  media  and  magnetic  pres¬ 
sure.  This  stage  is  characterized  by  inertial  axissymmetric  expan¬ 
sion  of  a  plasma  plume  and  almost  spherical  expansion  of  a  plasmoid 
moving  with  any  drift  velocity.  In  other  words  the  artificial  plas¬ 
ma  at  this  stage  does  not  "feel"  neither  magnetic  field  nor  envi¬ 
ronment  so  that  only  inner  electric  fields  control  the  plasma  flow. 

In  gas  releases  non-elastic  collisional  processes  including 
ionization  and  excitation  are  sufficient  at  the  first  stage  of  ex¬ 
pansion. 

The  second  stage  of  flow  expansion  corresponds  to  the  conditi¬ 
on  when  magnetic  pressure  already  exceeds  that  one  of  plasma  where¬ 
as  the  density  of  the  latter  is  still  sufficiently  greater  than  en¬ 
vironmental  density  so  that  the  movement  of  particles  conserves 
collective  character  and  plasma  expansion  across  the  magnetic  field 
is  inhibited.  As  a  result  the  plasma  plume  obtains  a  form  of  three 
dimensional  "petal"  and  the  plasmoid  becomes  "a  knitting  needle" 
stretched  along  the  magnetic  field. 

At  the  third  stage  the  movement  of  plasma  particles  becomes 
non-collective  and  is  completely  controlled  by  external  natural  fi¬ 
elds  and  interaction  with  environmental  particles.  This  region  may 
be  treated  as  a  "trace"  of  the  injection  which  used  "to  live"  for  a 
relatively  long  time.  This  is  a  weak  but  large  inhomogenity  usually 
interesting  for  geophysists  which  love  active  experiments  in  space. 
One  may  expect  that  this  region  is  also  responsible  for  return  flu¬ 
xes  of  injected  particles  to  the  spacecraft.  We  have  made  a  brief 
analyses  on  dynamics  of  the  third  stage  of  plasma  expansion  in  con¬ 
nection  with  interpretation  of  certain  our  previous  space  experi¬ 
ments.  This  approach  if  being  interesting  could  be  also  a  base  of 
future  cooperation. 


For  typical  injections  of  modest  intensity  (flow  rates  of  1019 
-  1021  particles  per  second  and  102°  -  1022  per  pulse)  the  charac¬ 
teristic  dimensions  of  the  first  stage  of  flow  ranges  up  to  ten  me¬ 
ters,  of  the  second  stage  -  from  a  few  tens  to  a  few  hundreds  me¬ 
ters  and  of  the  third  stage  -  from  a  few  hundreds  to  a  few  thousand 
meters. 

Our  experience  (and  not  only  ours)  suggests  that  attempt  to 
formulate  unitary  dynamic  problem  for  all  stages  of  APF  expansion 
results  in  large  non-visual  and  inconvenient  algorithms  of  enormous 
run  times.  Much  more  effective  way  is  formulation  and  solution  of 
separate  dynamical  problems  for  every  stage  of  expansion  of  every 
APF  type. 

Such  an  approach  results  in  relatively  compact  visual  soluti¬ 
ons  suitable  either  for  "intuitive"  analyses  or  numerical  calcula¬ 
tion. 

While  speaking  about  dynamical  and  mathematical  models  of  APF 
we  actually  separates  two  distinct  phases  in  the  procedure  of  dyna¬ 
mic  problem  formulation  and  solution. 

The  first  phase  encloses  evaluation  of  acting  forces  and  "dy¬ 
namical  range"  of  the  problem,  consequent  reduction  of  the  initial 
set  of  equations  to  a  final  dimensionless  form  convenient  for  solu¬ 
tion  . 

The  second  phase  results  in  an  analytical  solution  of  the 
problem  and/or  in  a  code  for  numerical  calculation.  Surely  it  is 
the  general  routine  of  dynamical  problem  solution.  With  this  sepa¬ 
ration  of  the  routine  in  two  phases  we’d  like  to  underline  that 
principally  there  is  or  may  be  a  few  ways  to  obtain  the  mathemati¬ 
cal  result  basing  on  the  same  or  similar  dynamical  model. 

In  the  ranges  of  our  interests  the  APF  mathematical  models  are 


mainly  used  for  description  and  study  of  radio  frequency  (RF)  ref¬ 
raction/scattering  effects.  The  base  of  scattering  calculations  is 
the  ordinary  method  of  the  eikonal  equation  numerical  solution  with 
dielectric  permeability  being  expressed  via  electron  density  dist¬ 
ribution.  Method  of  linearization  for  the  eikonal  initial  set  of 
characteristic  equations  now  under  development  allows  for  suffici¬ 
ent  reduction  of  run  time. 

For  every  APF  type  and  for  every  stage  of  expansion  we’d  like 
to  present  the  following  PC  codes: 

*spatial  intensity  distribution  of  scattered  spherical  wave 
(usable  example  gratias  for  evaluation  of  communication 
conditions  along  the  line  "board  -  the  Earth"); 

**spatial  intensity  distribution  of  scattered  plane  wave 
(the  same  for  the  line  "the  Earth  -  board"); 

***cross  section  of  backward  scattering 

(usable  for  evaluation  of  observation  conditions). 

Although  the  similar  approaches  are  used  in  development  of  all 
these  codes,  every  one  is  the  separate  independent  program  product 
since  it  includes  its  own  APF  model  and  utilizes  its  own  way  of  op¬ 
timization. 

Further  symbols  *,  **,  ***  are  used  to  denote  the  code  type. 


[ 


The  works  mentioned  above  may  be  displayed  with  the  scheme: 


*  ** 


This  scheme  also  shows  our  proposal  on  annual  distribution  of 
the  works  mentioned  above  with  account  for  the  contracts  In  force 
and  actual  volume  of  work. 

Finally  we  propose  the  following  annual  distribution  of  the 
works. 


Contracts  of  1994  in  force 


■25000$ 

workl-94 

SPC-94-4081 

Dynamical  and  mathematical  models  of 

plume  and  plasmoid  dense  core  ; 

20000$ 

work2-94 

SPC-94-4080 

Dynamical  and  mathematical  model  of 

ionized  jet. 

-  o  - 

•23000$ 

Proposed  contracts  of  1995 

work3-95  Dynamical  and  mathematical  model  of 

3-dimensional  plasma  petal; 

•23000$ 

work4-95 

Linearization  of  eikonal  equation 

and  code  *  for  plume  dense  core; 

Pronosed  contracts  of  1996 

work5-96  Codes  *  and  **  for  plume  dense  core; 

work6-96 

Codes  *,  **  and  *** 

for  3-dimensional  plasma  petal; 

(This  contracts  being  fulfilled  close  the  vertical  chains  on  plu¬ 
me  dense  core  and  3-dimensional  plasma  petal  the  scheme  presented 
above) . 


Proposed  contracts  of  1997 
work7-97  Dynamical  and  mathematical  model  of 
3-dimensional  knitting  needle; 
work8-97  Codes  *,  **  and  ***  for 

3-dimensional  knitting  needle; 
work9-97  Codes  *,  **  and  **  for  ionized  jet; 

(This  contracts  being  fulfilled  close  the  vertical  chains  on 
plasmoids  and  ionized  jets). 


At  the  bottom  line  we  propose  two  works  for  1995  -  one  on  dy¬ 
namics  and  one  on  scattering  ,  also  two  works  for  1996  -  both  on 
scattering  and  three  works  for  1997  -  one  on  dynamics  and  two  on 
scattering. 


